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Abstract 

We present the Monte Carlo event generator WINHAC for the charged-current 
Drell-Yan process, i.e. single VF-boson production with leptonic decays in hadronic 
collisions. It features multiphoton radiation within the Yennie-Frautschi-Suura 
exclusive exponentiation scheme with 0(a) electroweak corrections. For the initial- 
state QCD/QED parton shower and hadronisation it is interfaced directly to PYTHIA. 
It implements the options for proton-proton, proton-antipronton and any nucleus- 
nucleus collisions. Moreover, it includes the polarized VF-boson production and the 
neutral current Drell-Yan process, both at the Born level. Generation of weighted 
as well as unweighted (weight= 1) events is possible. WINHAC has been thoroughly 
tested numerically and cross-checked with independent Monte Carlo programs, such 
as HORACE and SANC; in the so-called tuned comparisons, an agreement at the 
sub-per-mil level has been reached. It has also been used as a basic tool for develop- 
ing and testing some new measurement methods of the Standard Model parameters 
at the LHC. 
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1 Introduction 



Motivations for investigating the Drell-Yan processes, in particular that involving the 
W-bosons, at the Tevatron and LHC can be found e.g. in Refs. [H [21 El [H El El [71 El El [10] . 

In theoretical descriptions of the Drell-Yan processes reaching a sufficiently high pre- 
cision for the Tevatron and LHC experiments requires including besides the QCD effects 
also the electroweak (EW) corrections, in particular the QED final-state radiation (FSR) 
P]. It was shown in Ref. [11] that including 0(a) EW corrections may be not sufficient 
because the higher-order FSR effects can shift Mw by ~ 10 MeV. 

Regions of a high transverse momentum of a charged lepton or a large W-boson 
transverse mass will be used at the LHC for various "new physics" searches. In such 
regions the O(o) EW corrections beyond the pure FSR can reach the size of 20-30% 
[T2l fl~3] . Given high statistics of VT-boson data expected at the LHC, these effects need 
to be included in order to have precise enough description of the Standard Model (SM) 
background. 

Drell-Yan processes will also be treated at the LHC as the so-called 'standard candle' 
processes (detector calibration, normalisation, etc.). Therefore, their precise theoretical 
description is very important. In order to be fully exploited by the experiments, such a 
description ought to be provided in the form of a Monte Carlo event generator. Several 
semi-analytical programs as well as Monte Carlo integrators and event generators for 
Drell-Yan processes including QCD and/or EW effects have been developed in recent 
years, see e.g. Refs. [HI [TO] • 

In this paper we briefly describe the Monte Carlo event generator WINHAC [H] 
dedicated to the charged-current Drell-Yan process. Section 2 is devoted to the QED 
FSR, while in Section 3 we describe inclusion of the 0(a) EW corrections. Finally, 
Section 4 contains summary and outlook. 



2 QED FSR 

Our first step in the development of WINHAC was to take into account effects of multipho- 
ton FSR [T5], based on the Yennie-Frautschi-Suura (YFS) exclusive exponentiation [TB] . 
The total cross section for a parton-level process 

gi(pi) + — W ± (Q) — l( gi ) + u(q„) + 7(h) + . . . + 7 (AvO, (n = 0, 1, . . .) (1) 
can be expressed as 
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The exponent in the above equation is the so-called YFS form factor in which all the in- 
frared (IR) singularities coming from virtual and real photons are re-summed (in a gauge- 
invariant way) to the infinite order and properly cancelled. The rest of the equation is 
the accordingly reorganized perturbative series, which in the current version of WINHAC 
is truncated at 0(a). In particular, the /3q and fti functions are the so-called non- 
IR YFS residuals corresponding to zero and one real hard-photon emission, respectively. 
This means that the former contains pure virtual corrections while the latter the pure 
real-photon ones. All details are given in Ref. [TB] . 

Based on the above formula an appropriate Monte Carlo (MC) algorithm was con- 
structed. We employed the importance sampling method by simplifying step-by-step the 
starting formula ([2]) until the simple enough probability density was obtained that could 
be generated with the use of basic MC methods. All these simplifications were compen- 
sated with corresponding MC weights. From these MC weights the full event weight was 
constructed. It can be used directly in weighted-event mode, e.g. for histogramming, or 
it can be applied in a rejection loop to construct weight= 1 events in unweighted-event 
mode. In the weighted-event mode the program calculates for each event a large array 
of the MC weights corresponding to various effects, corrections, etc. The parton level 
cross section is convoluted with a product of two parton distribution functions (PDFs) 
corresponding to two colliding beams: proton-proton, proton-antiproton or any nucleus- 
nucleus. Values of Bjorken-x's of a quark and an antiquark participating in the hard 
process are generated from a 2-dimensional distribution, given by the product of two 
PDFs convoluted with a simplified ("crude") parton-level cross section, with the help of 
the adaptive cellular MC sampler FOAM [17J. Then, the quark and antiquark flavours 
are randomly chosen according to their relative contributions to the cross section. 

After constructing the Monte Carlo event generator WINHAC we performed several 
numerical tests, both internal and external. The important external test was the com- 
prehensive comparison with the independent MC program HORACE [11] which includes 
higher-order FSR corrections through a QED parton shower algorithm. We found the 
agreement between the two programs at the per-mil level, see Ref. [TS] for details. This 
is very important since the two programs are based on very different approaches to the 
QED corrections. It suggests that major parts of the higher-order FSR corrections are 
included in both programs. 

3 EW corrections 

Strictly speaking, for the charged-current electroweak processes it is not possible to split 
the EW corrections into the pure QED and pure weak ones in a gauge-invariant way. 
However, one can find gauge-invariant subsets of QED corrections and separate them 
from the rest. Such a separation is in general not unique, as one can identify different 
subsets of gauge invariant QED corrections. In our definition of the QED FSR corrections 
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described the previous section we used the prescription of Marciano and Sirlin (M&S) 
given in Ref. [19]. Then we included the rest of the electroweak corrections for the FSR 
in the pole approximation (PA). These corrections have been provided to us in the form 
of the Fortran-code module by the SANC group [20]. Numerically they are small - at the 
level of 0.2-0.3% - and almost constant, so that they are not important for the ly-mass 
measurements. 

The main problem that appears in including the full 0(a) EW corrections (beyond 
PA) to the charged- current Drell-Yan process concerns the QED initial-state radiation 
(ISR), i.e. photon radiation from incoming quarks. In EW calculations the collinear ISR 
singularities are usually regularized with finite quark masses m q . As a results, terms with 
ln(s/mg) appear in theoretical predictions. However, m q in this context is not a well- 
defined quantity, mainly because of the QCD effects. In fact, these singularities can be 
absorbed into the PDFs in a similar way as the ones in QCD, which is usually done in 
the MS or DIS scheme, see e.g. Refs. [121 [131 [20] . QCD parton-shower generators, such as 
PYTHIA or HERWIG, usually generate also the QED ISR, however they do it in neither 
of the above schemes. So, combining such MS or DIS-subtracted EW corrections with 
parton shower generators would result in some mismatch. Besides, the above schemes are 
not well suited for fully exclusive Monte Carlo generators, as they are defined for partly 
integrated quantities. 

Therefore, we take a different approach to the QED ISR. We identify the QED gauge- 
invariant subsets of the ISR corrections and subtract them from the total 0(a) EW 
correction, leaving responsibility for their computation to parton-shower generators. In 
our opinion, it is better if the parton-shower generator treats both the QCD and QED 
ISR since the QED radiation is strongly affected by the QCD one. Of course, the parton- 
shower programs usually include the QED ISR in some approximation. However, the 
QED ISR corrections are small and flat for the main ly-boson observables [1] . Moreover, 
at the LHC the W-boson distributions will have to be normalized to the corresponding 
Z-boson ones in order to avoid some large uncertainties, e.g. coming from the PDFs 
parametrisations, ly-boson p^-spectrum predictions, lepton energy scale determination, 
etc. In such ratios most of the QED ISR uncertainties will also cancel. 

In WINHAC we have implemented three schemes for ISR subtraction: the one simi- 
lar to the M&S prescription for the FSR, the YFS-motivated scheme and the one from 
Ref. [21J. In fact, in the latter two schemes the ISR virtual plus soft real-photon correc- 
tions consist of the YFS form factor plus the terms ~ (l/2)Q^[ln(s/mf) — 1], however 
our own formulae differ from those of Ref. [21] by some constant terms. Similar schemes 
can also be applied to the FSR and initial-final interferences - resulting in a separation of 
the 0(a) EW corrections into the 'QED-like' corrections and the 'weak-like' corrections. 
Such a splitting is, of course, not unique and it differs in various schemes by some constant 
terms. Our numerical tests show that the YFS-motivated scheme for the ISR is very close 
to the MS scheme - the relative difference is < 0.05%. 

Similarly as in the case of the FSR, for the full O(a) EW virtual corrections we use 
the Fortran-code modules from the SANC library [20] . Both virtual and real-photon 
corrections are implemented on top of the formula ([2]) through appropriate MC weights. 
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In order to test their implementation, we performed numerical comparisons at O(o) of 
WINHAC and the SANC MC integrator based on VEGAS. We have found a very good 
agreement - at the sub-per-mil level - between two programs for the total cross sections 
and various distributions, see Ref. [22] for more details. 



4 Summary and outlook 

We have briefly presented the Monte Carlo event generator WINHAC for precision theo- 
retical prediction of the charged-current Drell-Yan process, i.e. single iy-boson produc- 
tion with leptonic decays in hadronic collisions. It includes multiphoton radiation in the 
YFS exclusive exponentiation scheme with the ISR-subtracted 0(a) EW corrections. For 
QCD/QED initial-state parton shower as well as hadronisation, it is directly interfaced 
with PYTHIA 6.4. It has been tested numerically and cross-checked to the high preci- 
sion with independent programs, in particular HORACE and SANC. It has been used 
for developing new strategies of measuring the SM parameters (51 [TJ El [9] and testing the 
mechanism of the electroweak symmetry breaking [3j at the LHC. 

The current version, 1 . 30, of WINHAC is written in F0RTRAN77, however work is 
now in progress on re- writing it in C++. A similar program, called ZINHAC [23], for the 
neutral-current Drell-Yan process is being completed now - in C++. 
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